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Abstract. We present a dynamical model for the formation and evolution of a massive disk galaxy, within a growing dark halo 
whose mass evolves according to cosmological simulations of structure formation. The galactic evolution is simulated with a 
new three-dimensional chemo-dynamical code, including dark matter, stars and a multi-phase ISM. The simulations start at 
redshift z = 4.85 with a small dark halo in a ACDM universe and we follow the evolution until the present epoch. The energy 
release by massive stars and supernovae prevents a rapid collapse of the baryonic matter and delays the maximum star formation 
until redshift z ~ 1. The metal enrichment history in this model is broadly consistent with the evolution of [Zn/H] in damped 
Lya systems. The galaxy forms radially from inside-out and vertically from halo to disk. As a function of metallicity, we have 
described a sequence of populations, reminiscent of the extreme halo, inner halo, metal-poor thick disk, thick disk, thin disk 
and inner bulge in the Milky Way. The first galactic component that forms is the halo, followed by the bulge, the disk-halo 
transition region, and the disk. At redshift z ~ 1, a bar begins to form which later turns into a triaxial bulge. Despite the still 
idealized model, the final galaxy resembles present-day disk galaxies in many aspects. The bulge in the model consists of at 
least two stellar subpopulations, an early collapse population and a population that formed later in the bar. The initial metallicity 
gradients in the disk are later smoothed out by large scale gas motions induced by the bar. There is a pronounced deficiency of 
low-metallicity disk stars due to pre-enrichment of the disk ISM with metal-rich gas from the bulge and inner disk ("G-dwarf 
problem"). The mean rotation and the distribution of orbital eccentricities for all stars as a function of metallicity are not very 
different from those observed in the solar neighbourhood, showing that early homogeneous collapse models are oversimplified. 
The approach presented here provides a detailed description of the formation and evolution of an isolated disk galaxy in a 
ACDM universe, yielding new information about the kinematical and chemical history of the stars and the interstellar medium, 
but also about the evolution of the luminosity, the colours and the morphology of disk galaxies with redshift. 
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1. Introduction 

During the last decade, significant progress has been made 
in understanding cosmic structure formation and galac- 
tic evolution. With high-resolution cosmological Simula- 



tions, the formation of dark halos in different cosmolo- 


gies has been studied in detail (e.g. 


Navarro et ak, 


1996 


Moore et all 1 19981 ICoKn et all 12000: 


Yoshida et all 


2000 


Avila-Reese et all l200ll Ijenkins et all 


2001; Klvmnetal. 


2001) to determine the large-scale mass distribution, the 



halo merging histories, the structural parameters of the 
dark halos, and fin ally the formation of g a laxies in- 
side t h ese dark halos (j Steinmetz & Mulle n 1 19951 iKav et all 
200ri iNavarro & Steinmet A l2000bl iMosconi et all 1200 U 
Pearce et alll200ll) . 

In many respects, these simulations are in agreement 
with observations, but there remain a number of diffic ulties. 
Firstly, the universal cuspy halo density profiles dNavarro et all 
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Il997t iMoore et all Il999t Fukushige & Makino, 2001 ), while 
in agreement with observations of galaxy clusters, are in 
apparent disagreement with the flat dark matter density 
distributions infer red from observations in the centres of 
galaxies dSalucci & BurkertLl200filBlais-Ouellette et allhoOlt 
iRorriello & SaluccilkOOlUde Rlok et alll200ll) . Secondly, the 
specific angular momenta and scale lengths of the disk galax- 
ies in the cosmological simulations are too small compared to 
real galaxies (Navarro & Steinmetz, 2000b). This may be par- 
tially a numerical problem of SPH simulations, but is more 
likely due to overly efficient cooling and strong angular mo- 
mentum trans port from the baryons to the dark halos du ring 
merge r events dNavarro & BenzUl99ll).lThacker & Couch man 
(2001) and also Som me^Larsenet alJdl999i) showed the angu- 
lar momentum problem in CDM simulations may be resolved 
by the inclusion of feedback processes. Thirdly, the cold dark 
matter (CDM) simulations predict a large number of satellite 
galaxies in the neighbourhood of gian t galaxies like the M ilky 
Way Galaxy, but this is not observed dMoore et allll999T) . 
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The cosmological simulations are used to describe the 
large-scale evolution of the dark matter distribution from the 
primordial density fluctuations to the time when the dark 
halos form. It is inevitable in such simulations that struc- 
tures on galactic and sub-galactic scales are not well-resolved. 
Furthermore, it is difficult to incorporate the detailed physics 
of the baryonic component, that is of the multi-phase interstel- 
lar medium (ISM), and the feedback processes between stars 
and ISM. These are serious drawbacks, because feedback from 
stars can prevent a proto-galactic cloud from rapid collapse an d 
it can trigger large scale gas motions (Samla nd et all 11997). 
Both alter the galactic formation process significantly. 

Complementary ways to circumvent at least some of 
these problems are (i) to combine th e large-scale simulations 
with semi-analytical g al axy m odels (Guiderdoni et al., 1998; 
Somerville & PrimacM 1 19991 iBoissier & Prantzos. 2000; 



Cole et all bOOCt iKauffmann & Hae hneltt l2000l) . or (ii) to 



simulate the formation of only single galaxies or small 
galactic groups using either cosmological initial conditions 
or a simplified, but cos mologically motivated model of the 
dark matter background ( Navarro & Steinmetz, 1997; Berczik, 
1999; Hultman & Pharasvnl 1 19991: ISommer-Lars en et all 
19991 iBekki & Chibal I200U llTiacker & Couchmanl 120*011 
Williams & NelsonL 1200 lb .While the semi-analytical models 



are advantageous for investigating the global properties of 
galaxy samples, the small-scale dynamical models provide 
information about the detailed structure of galaxies, the 
kinematics of the stellar populations and the ISM, and the 
star formation hist ories. The dynam i cal models are called 

j — 1 r — — 1 1 1 r~ — 

chemo-dynamical (Theis et al., 1992; Samland et al., 1997; 
Bek ki & Shiovl Il998t iBerczik L Il999t IWilliams fe Nelsonl 
2001), if they include different stellar populations, a multi- 
phase ISM, and an interaction network that describes the mass, 
momentum and energy transfer between these components. 

The aim of the present work is to investigate how a large 
disk galaxy forms inside a growing dark halo in a realistic 
ACDM cosmogony. We present the results of a new three- 
dimensional model, including dark matter, stars, and a multi- 
phase ISM. The mass infall at the boundaries of the simulated 
volume is tak en from cosmological simulations (see VIRGO- 
GIF-project; Kauff mann et allfl999l) . We use a total dark mat- 
ter mass of 1.5xl0 12 M©, a total baryonic mass of 3x 10 11 Mq, 
a spin parameter A = 0.05, and an angular momentum profil e 
similar to the u niversal profile found by lBullock e t all J200ll) . 
Different from IBekki & Shioval Jl998t) : iBerczikl Jl999h and 
IWilliams & Nelson! J200 lh . we use an Eulerian grid-code for 
the ISM, which we believe is better suited to describe the 
multi-phase structure of the ISM and the feedback processes, 
especially in low-density regions. The final galaxy model pro- 
vides densities, velocities, velocity dispersions, chemical abun- 
dances, and ages for 614 500 stellar particles, and temperatures, 
chemical abundances, densities, velocities, and pressures of the 
ISM. The goal of the project is a self-consistent model for the 
formation and evolution of large disk galaxies, which can be 
used to understand observations of young galaxies as well as 
data on stellar populations in the Milky Way. 

The paper is organized as follows. Section |2] summarizes 
some observational material relating to the formation and evo- 



lution of galaxies. Section [3] describes our multi-phase ISM 
model and the underlying dark matter halo model. In Section|4] 
we outline the formation history of the model disk galaxy. The 
properties of the galactic halo, bulge and disk components are 
discussed in Section [5] A summary and a short outlook con- 
clude the paper (Section|6j. 

2. Observations of the galaxy formation process 

Galaxy evolution models can give a simplified description of 
the evolution of real galaxies and can make predictions that can 
be tested by observations. The high redshift observations show 
that the first proto-galactic structures form at redshift z > 3 (e. 
Clements & Coiichlll994lDev et allfl998tlSteidel et all fTq 



Manning et alll2QO0l) . Surveys of high redshift galaxies reveal 
a wide range of galactic morphologies with considerable sub- 
structure and dumpiness JPentericci et all l200ll) . At redshift 
z ~ 1-2, the formation of galactic components (sphe roids , 
halos , bulges and disks) s eems to take place JLillv et all 
1998; iKaiisawa & Yamadal l200ll) . However, the spatial res- 



olution of these high and intermediate redshift surveys (e.g. 


Abraham et al. . |l999|: 


van den Bergh et all l2000l lEllis et all 


EoOlUMenanteau et al. 


l200ll) is still too low to measure more 



than global galactic properties (e.g. asymmetry and concentra- 
tion parameters). There are indications th at most of the massive 
galax ies form before or around z — 1 JBrinchmann & E llis. 
200o]) , and especially th at the gia nt ellipticals form at reds hifts 
z > 1 JBest et allll998l). (but see Ivan Dokkum et ail fl999). but 
not before z = 2.4 dPentericci et all 1200 ll) . There is a general 
trend of early type galaxies being older and h aving higher pea k 
star formation rates than late type galaxies (Sandage, 1986). 
These general star formation histories seem to be overlapped 
with star formation bursts, triggered by infalling material or in- 
teractions with other galaxies. 

Since we can observe galactic sub-structures down to the 
stellar scale only in our Galaxy, the Milky Way observations are 
a cornerstone for the understanding of the formation and evo- 
lution of all spiral galaxies. From studies of Galactic stars we 
know that the halo is old and that the disk and bulge contain a 
mixture of old and young stars. This raises the question in what 
sequence the galactic components halo, bulge and disk formed, 
and what was the star formation history. Determinations of 
the star formation history of the Milky Way and other nearby 
galaxies were attempted bvlBell et all ( l2000l) : lHernandez et all 
(2000]) and iRocha-Pinto et all J2000I) . based on star counts, 
broad band colours, or chemical compositions of stars. 

The chemical and kinematical data from galactic halo stars 
have been one of the starting points for dif ferent scena r ios for 
the formation and evolution of galaxies. Egge n et all d!962l) 
were the first to propose a monolithic collapse galactic for- 
mation model in which the Milky Way forms rapidly out of 
a collapsing proto-galactic cloud. This was based on data for 
galactic halo stars which are now known t o have been plagued 
by selection effects (Chiba & Beers, 2000). Later, based on ob- 
servations of halo globular clusters. ISearle & Zinnl (1 1978) pro- 
posed a new model in which the Milky Way forms out of merg- 
ing lumps and fragments which all have their individual star 
formation histories. The contrasting nature of these two models 
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and the proposed observable differences (e.g. existence or non- 
existence of metallicity gradients, the expected age spread for 
halo objects) seemed to offer a simple way to deduce the for- 
mation history of disk galaxies from Milky Way observations. 
However, new data especially from proper motion catalogues, 
lead to the conclusion that at least the halo of the Milky Way 
formed neither by a monolithic collapse nor by a pure merging 
proce ss JSandagetfl990tlUnavane et allll996tlchiba & Beerl 
2000). In the modern hybrid scenarios, a part of the Galactic 
halo formed during a dissipative process (collapse) and an- 
othe r part show s signatures of merging events. For example, 
IChiba & Beers! (EoOCO proposed a scenario in which "the outer 
halo is made up from merging and/or accretion of sub-galactic 
objects, such as dwarf-type satellite galaxies, whereas the inner 
part of the halo has undergone a dissipative contraction on rel- 
at ively short time sca les". This is consistent with the findings 
of iHelmi etail dl999h . who found that 10% of the metal-poor 
stars in the outer halo of the Milky Way are aligned in a sin- 
gle coherent structure, which they identify as the remnant of a 
disrupted dwarf galaxy. These hybrid scenarios are, in the gen- 
eral outline, also consis tent with the f ormation of galaxies in a 
hierarchical universe (Chiba & Beersl EoOO). 



3. The model 

The formation and evolution of a galaxy is influenced by the 
cosmological initial conditions and the environment, but also 
by internal feedback processes. Heating by supernovae (SNe), 
dissipation, radiative cooling, as well as in- and outflows influ- 
ence the evolution both locally and globally. Our new three- 
dimensional galactic model takes into account the dynam- 
ics of stars and the different phases of the ISM, as well as 
the processes ("chemistry") which connect the ISM and the 
stars. These baryonic components are embedded in a collision- 
free dark matter halo taken from simulations of cosmological 
structure formation. The description of the ISM is based on 
that used in the two-dimensional chemo-d ynamical models of 
ISamland & Henslerl ill 99(1) : ISamland et alJlll997l) . The follow- 
ing subsections contain a brief outline of our model. 



3.1. The dark halo 



Cosmological simulations (e.g. [ 


^avarro et al l 


Moore et all 


1999; 




1999 


) show that 



1995; 
l a hi- 
erarchical universe galaxies form out of a few large fragments 
with an additional slow accretion of mass. Thus each galaxy 
has its own formation history and characteristic shape, 
depending on the masses, angular momenta and merging 
histories of the fragments. Yet the formation of the CDM 



IWechsler et alJ J20021 and 


vanden Boschl J2OO2I) calculated 


the mass accretion history, 


Bullock et a 


. (2001) the angular 


momentum distribution. Navarro et alJ i 


1997) 


and others the 


density profiles and Del PodoIoI 


2001 


and IWechsler et al. 




(2002) the evolution of the concentration parameter. 



time [Gyr] 



Fig. 1. The dots shows the average mass from 96 different halo 
merging trees, which all result in a 1.8 x 10 12 Mq halo. The 
smooth line shows the approximation which is used in our 
numerical model. The most extreme halo merging trees and 
that which shows the smallest deviations from the average are 
shown in the smalle r upper pane ls. Data from the VIRGO-GIF- 
project JKauffmann et all 1 19991) 



3.1 .1 . The dark matter halo formation history 

In the hierarchical formation scenario, a dark halo may form 
early and rapidly or grow slowly with time. For the model de- 
scribed in this paper we take an average halo formation his- 
tory, because at present we cannot simulate the evolution of 
a galaxy for a large number of different halo formation histo- 
ries. We choose a ACDM universe with A = 0.7, Q = 0.3, a 
Hubble constant ho = 0.7 and a baryonic-to-dark matter ratio of 
1 :5. From the cosmologica l N-body s imulations of the VIRGO- 
GIF-project JKauffmann et all |l999), we extracted the merg- 
ing trees of 96 halos, all with a final mass of 1.8 x 10 12 Mq. 
Figure ^ shows the average halo mass growth, which is in 
agreement with the un iver sal mass accretio n histories found by 
Ivan den Boschl J2002I) and lWechsler et alJ J2002I) . The individ- 
ual merging histories can be quite different; two extreme cases 
with early and late mass growth are shown in the small panels 
at the top left and right of Fig. ^ 

We assume that the gro wing dark halo has t he universal 
density profile as proposed bv lNavarro et alJ Jl997l) . and that its 
mass is accreted in spherical shells. Given the current halo mass 
and the evolution of the criti cal density, we can determine the 
outer halo radius r2oo (see iMo et all fl998). A free parameter 
that remains is the mass concentration c, which in the present 
simulation is assumed to be c = 9.2. At redshift z — 0, this leads 
to a circular velocity of v c j rc = 179 km s at the /-200 = 250 kpc 
radius. Inside a galactocentric radius of 8.5 kpc, the final mass 
of dark matter is 4 x 10 10 M0, comparable to the upper limit of 
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the d ark matt er content insid e the solar circle (5.2 x 10 1() Mq; 
iNavarro & Steinmetz , 2000a). 

Treating the growing dark matter halo as a spherical back- 
ground potential that does not respond to the baryonic compo- 
nent minimizes gravitational torques during the collapse and 
ther eby circumvents the well-know n angular momentum prob- 
lem (tNavarro & Steinmetz, 2000b). The accretion in this model 
is smoother as compared to hierarchical models, where a frac- 
tion of the mass accretion takes place through merging dark 
matter fragments, unless the gas in these hierarchical fragments 
is dispersed by early star formation and feedback. 



1000 



100 z 




20 30 
r [kpc] 



3.1 .2. The angular momentum of the dark matter 

During the hierarchical clustering, tidal torques induce shear 
flows which lead to rotating proto-galactic systems. N-body 
simulations show that, independent of the initial perturbation 
spectrum and the cosmological model, the spin parameter of 
the dark halos is in the range A = L\E\ {) - 5 G ' 1 M~ 2 - 5 _ 0.02-0.11 
with an average value of A = 0.0 4 — 0.05 jBarnes & Efstathiou} 
19871 ISteinmetz & Bartelmanrl 1 199.4 ICole &Lacevi 1 199(4 
van den Boschlll998UGardnerll200llK 

Bull ock et alJ ( 1200 II) studied the angular momentum pro- 
files of dark halos in a ACDM cosmology. They found that the 
spatial distribution of the angular momentum in most of the 
halos show a certain degree of cylindrical symmetry, but that 
a spherically symmetric angular momentum distribution with 
a power-law j(r) oc r 11±0 - 3 i s also a good approximation. We 
consider the following two rotation fields for the halo: 



halo _ 
v rot ~ V (T 



^x 2 + y 2 + z 2 + do 



= v'o 



\d\ + d 



. .halo 



v 



V-* 2 + y 2 
^x 2 +y 2 + d 



= vo- 



■ do 



(1) 



(2) 



vo and do are constants which determine the total amount of an- 
gular momentum and the rotation curve near the rotation axis, 
respectively. 

Assuming that initially the mass and velocity distribu- 
tions of the baryonic matter coincide with those of the dark 
matter, and that baryonic mass and specific angular m omen- 
tum is cons erved during the collapse to a disk (Mestel ll963t 
iFall & Efstathiou, 1980), we can estimate the disk surface den- 
sity profiles corresponding to these rotation fields. Figure |2] 
shows the results, calculated for c ollapse in the pur e NFW 
potential. As already pointed out by Fall & Efstathiou ( 1980), 
disks with exponential surface density profiles are found when 
the halo has a rising rotation curve at small radii and a nearly 
constant rotation velocity in the outer parts. In the present case, 
the cylindrical halo rotation model (eqn. |5Jl also produces an 
exponential disk, at least in the inner 20 kpc. The spherical ro- 
tation model (eqnQ} instead leads to higher surface densities 
near the galactic centre which favours the formation of a bulge. 
In the following we use the spherical rotation model with a spin 
parameter A = 0.05, to fix the initial angular momentum distri- 
bution of the baryonic matter. 



Fig. 2. Surface density of the galactic disk if the specific angu- 
lar momentum is conserved during the proto-galactic collapse 
in our NFW halo model. The cylindrical rotation leads to an 
exponential disk in the inner 20 kpc (dashed curve), while the 
spherical rotation produces an exponential disk with an addi- 
tional mass at the centre (full line). The thin dash-dotted line 
shows an exponential disk with a scale length of 8 kpc. 



3.2. The Stars 

The stellar motions have a considerable influence on the galac- 
tic evolution, since the time-delayed stellar mass loss and en- 
ergy release connect different regions and also different galac- 
tic evolution phases. The approximately 10 11 stars in a large 
galaxy form a collision-less system. Each star moves on its own 
orbit, which is determined by the conditions under which the 
star formed and the evolution of the g alactic gravitational po- 
tentia l. We use a particle-mesh method JHocknev & Eastwood! 
1988), in which the passing of the stellar particles to the mesh 
and also the interpolation of the gravitational forces at the parti- 
cle positions is done with a cloud-in-cell method. This method 
is fast and second order accurate, but it is limited by the spatial 
resolution of the underlying grid, which at present is 370 pc. 

Each stellar particle represents an ensemble of stars with 
a stellar mass function. In the present model, we assume a 
Salpeter initial mass function (IMF) with lower and upper mass 
limits of 0.8 M0 and 50 Mq, and with an additional lock-up 
mass fraction of 60%. The lock-up mass equals the mass of 
all stars below 0.8 Mq. During a simulation the star forma- 
tion rate is calculated on the spatial grid, and in timesteps of 
lOMyr representative samples of 500 new stellar particles are 
created. This yields 614500 particles at the end of the simula- 
tion. Depending on the star formation rate a single particle rep- 
resents 10 4 to 10 7 stars. The initial chemical composition and 
velocity of a stellar particle is derived from the chemical com- 
position, velocity and velocity dispersion of the parent molec- 
ular cloud medium. 

We assume that the evolution of the stars is determined only 
by their initial mass and metallicity, neglecting effects of stel- 
lar rotatio n or star-star interactions. From the stellar evol ution 
models oflMaeder & Mevnet Nl989h :ISchaller et al. ( 1992) and 
ISchaerer et alJ Jl993l) we get the lifetimes, energy releases and 
mass losses on the main sequence. We can subdivide stars into 
three main classes. 
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First, the low mass stars which stay on the main sequence 
during the whole galactic evolution. These stars become notice- 
able only by their gravitational forces and because they lock-up 
a major fraction of the galactic mass. The total lock-up mass 
has influence on the galactic gas content and therefore on the 
sequence of stellar generations. It is important in smaller galax- 
ies where SN-driven winds reduce the mass in gas, but not the 
mass locked up in low mass stars. 

The second class of stars are intermediate in mass. When 
they leave the main sequence, they loose a substantial fraction 
of the initial mass (asymptotic giant branch and planetary neb- 
ula phase); however, the energy release of these stars is small 
compared to other energy sources. 

The third class of stars are the massive stars, typically with 
masses in excess of 8-10 M©. These stars are a major source of 
both energy and mass return in a galaxy. The UV-radiation and 
the kinetic energy release during the final explosion (SN type 
II) ionize and accelerate the ISM. This can trigger large-scale 
gas flows and star formation which both influence the galactic 
evolution. 

A rare but important galactic event is the explosion of SNe 
of type la. The most probable progenitor candidates are close 
binary stars consisting of a white dwarf and a main sequence 
star. We include th is process, since i t is a major source of iron 
peak elements dNomotoet allll984 . 

3.3. The Interstellar Medium 

The ISM in galaxies exists at different densities and tempera- 
tures, and it shows phase transitions on short timescales com- 
pared to the age of a galaxy. In the present galactic model 
we use a simplified ISM description based on the three-phase 
model of lMcKee~& OstrikeJ Jl977l) . The three phases are a 
cold, a warm and a hot phase, with typical temperatures of 
80 K, 8000 K and 10 6 K, respectively. The cold phase is found 
in the dense cores of molecular clouds. These cores are em- 
bedded in warm neutral or ionized gas envelopes. Outside the 
clouds the space is filled with hot, dilute gas. The three-phase 
model is a simplification, but it is realistic enough to use it in 
a global galactic model. A critical discussion of the w eak and 
strong points of the three-phase model can be found in McKee 
(1990). 

Essentially, the temperatures and densities of the ISM 
phases determine most of the physical processes in the ISM. 
However, the geometry of the clouds is also important for pro- 
cesses that take place in the phase transition regions (e.g. evap- 
oration and condensation) or that depend on the sizes of the 
clouds explicitly (e.g. cloud-cloud collisions and star forma- 
tion). In our model, we assume spherical clouds which follow 
the mass-radius relation of lElme greenl (H 9 8 9 ) 



i-»0)vs 



(3) 



P /k 

P4 — 1() 4^ cm j is the pressure in the ambient ISM. The mean 
cloud mass is 10 Mrr>; this is derived from the cloud mass 
spectrum of Da me et al I d 19861) for star forming giant molec- 
ular clouds in the mass range of 10 4 Mq to 10 6 Mq. 



In the present dynamical description, we use only a two- 
phase model for the ISM, containing hot gas with embed- 
ded (cold+warm) clouds. The motions of both components 
are described by the time-dependent hydrodynamical equa- 
tions (mass, momentum and energy conservation). Simila r as in 
the tw o-dimensional chemo-dynamical models (Sam land et all 
1 19971) we use a fractional step method to split the problem 
into source, sink and transport steps. The problem of the three- 
dimensional transport is broken down into a number of one- 
dimensional problems for the densities, metallicities, momenta 
and energies (di mensional sp litting). Using the van Leer ad- 
vection scheme (van Leer, 1977) combined with the consistent 
adve ction method INorman e t al.. 1980) and the Strang split- 
ting (Strang, 1968), the advection is second order accurate. In 
the Strang splitting the three-dimensional transport is split into 
a sequence of five one-dimensional transport steps: 



Af At At At 

AtA = (-A x )(-A y )(AtA z )(-A y )(-A x ) 



(4) 



where At is the time step length, A x , A y , A z are the differ- 
ential operators which describe the transport in x, y, z direc- 
tion, and A is the operator for the three-dimensional problem. 
Approximately the same accuracy can be achieved by a simpler 
splitting scheme which uses an alternating sequence of only 
three one-dimensional transport steps: 

(AtA)(AtA) = [(AtA x )(AtA y )(AtA z )] 

[(AtA z )(AtA y )(AtA x )] (5) 

For the ~ 10 6 transport steps in our galactic evolution simula- 
tions, both splitting schemes lead to identical results. However, 
the second method is faster by a factor of 1.6. 

3.4. The interaction network 

The interactions between stars and ISM are very important for 
the galactic evolution. They determine the evolution timescale 
of a galaxy to a large extent. In the following subsections we 
give a short overview about the processes which we take into 
account, and discuss the influence of the free parameters in the 
parametrization of these processes. 

3.4.1. Star formation 

Star formation is obviously one of the most important processes 
during the galac tic evolution, bu t it is one of the least under- 
stood. Early on, Schmidt ( 1959) found that the star formation 
rate in the solar neighbourhood varies with the squ are of the 
surfac e gas density 2 gas of the ISM. More recently, Ken nicutll 
( 1998) found a star formation rate for disk and starburst galax- 
ies which is proportional to X ga 4±0 - 15 , with a sharp decline 
below a critical surface density threshold (« 3-13 M pc~ 2 ). 
Remarkably, this simple law describes star formation in galax- 
ies for which the average gas consumption time can vary from 
3 ■ 10 8 yr (starbursts) to 2 • 10 9 yr (normal disks). 

Motivated by these findings we use the following star for- 
mation law: 

Pcld _ Pcld 



Psf = 



C s f?d 



(6) 



yn 
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Here p c \& is the average density of the cloudy medium, t S f the 
star formation timescale, and fdy n the dynamical timescale for 
the average star forming cloud mass. For the cloud model de- 
scribed in Section 1X31 the dynamical timescale of the average 
star forming cloud is 



t&y 



3n 



32Gp c id 



16.2 P 



-3/8 



(7) 



Molecular clouds do not collapse and form stars on a single 
dynamical timescale, because magnetic fields or turbulent mo- 
tions stabilize the clouds. To account for these processes, we 
include the factor c s f. If we take into account that a part of 
the clouds is only confined by pressure and that the timescales 
for ambipolar diffusion and dissipation of turbulent energy are 
much longer than the dynamical timescale, we expect c s f » 1 . 
c s f can be estimated fro m the observe d number of OB stars in 
the solar neighbourhood jReedll200lh . Assuming that the stel- 
lar mass function follows a Salpeter law with lower and upper 
mass limits of 0.1 Mq and 50 Mq, we derive a local star for- 
mation rate of 1.8 x 10~ 5 M pc~ 3 Myr~'. From the local ISM 
density and pressure (p « 0.035 M Q pc~ 3 ss 1cm -3 , P4 ss 1), 
we find c s f = 120, which means that the on average the star 
formation timescale exceeds the dynamical timescale by factor 
of 120. The final star formation law can then be written as 



Psf = ^Pcld^ 



3/8 



(8) 



with a star formation (gas depletion) timescale t s f ~ l/i] s f = 
1.9 Gyr in the local Galactic disk. The star formation rate is 
proportional to p l cl ^ 5 if the cloudy medium is in pressure equi- 
librium (P4 cc pdd). This is the quiescent star formation mode 
in the model. During the collapse or if large fragments merge, 
the pressure can increase and the star formation can be more 
efficient. 

Figure [3] shows the resulting star formation rate as a func- 
tion of the surface density of gas, in the forming disk galaxy 
model described in Section^] Both the average star formation 
rates and the star formation rates in the strongest star formation 

I . L 

regions are slightly higher than, but consistent with Kennicutt's 
1998 data and mean relation, over a range of surface densities 
extending from normal disk galaxies into the starburst galaxy 
region. Had we calibrated the factor 77^ on Kennicutt's diagram, 
the resulting value would have been about a factor of 2 smaller 
than that based on the solar neighbourhood OB stars. The typ- 
ical scatter in Kennicutt's plot corresponds to a factor of 2.2 
dispersion in r/ s f . 

3.4.2. Massive stars and supernovae type II 

Massive stars and SNe influence the galactic evolution in sev- 
eral ways. They heat the ISM, produce shells and filaments and 
stir up the cloudy medium. In addition, they are the most impor- 
tant source of heavy elements and they determine the chemical 
evolution of a galaxy. The lifetimes of the massive stars range 
from 3 to 30 Myr before they finally explode as SNe II. Only 
a stellar remnant of « 2 M0 remains. The rest of the mass is 
blown away by stellar winds or expelled during the explosion. 
On average a massive star on the main sequence emits about 
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Fig. 3. Star formation rate in the forming disk galaxy model 
described in Section [4] versus the projected surface density of 
the cloudy medium in a 20kpc sphere. Open numbered cir- 
cles are for average star formation rates and surface densities 
within 20 kpc; the number denotes the galactic age in Gyr. The 
filled numbered circles are for the most prominent star forma- 
tion region at the same galact ic ages. The black dots and black 
squares denote the data from lKennicuttl i 199 8) for disk galax- 
ies and starburst galaxies. The dashed line is Kennicutt's mean 
relation. 



6 x 10 51 erg in photons dSchaerer & Vaccalll998l) which ionize 
the surrounding ISM and regulate the star formation. During 
the SN explosion the star releases another 10 51 erg, creating a 
fast expanding bubble. This influences the star formation pro- 
cess, since it decreases the gas density in the surroundings of 
the SN and it can trigger star formation in other regions. 

In the model we can use only a simple description of these 
processes. 

(i) We assume that the mass return of a massive star takes 
place instantaneously at the end of its evolution, when the SNe 
explodes. This is a good approximation, because the time be- 
tween the stellar mass loss by a wind and the ejection of the SN 
shell is short compared to galactic timescales. 

(ii) We assume that the UV photons emitted by the massive 
stars are absorbed and reradiated by the dense cloud medium. 
Since the cooling time of this gas is short compared to the life- 
time of the massive stars (the typical timescale which we re- 
solve in our models), we do not simulate this heating and cool- 
ing process, but assume that the clouds are in thermal equilib- 
rium. 
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(iii) The SN explosion energy (10 51 erg) is released in a 
short time and produces a bubble of hot gas and an expanding 
shell. The shell interacts with the clouds and by this increases 
the velocity dispersion <x cld of the cloudy medium. The SN en- 
ergy is released locally, and we assume that rji cm = 95% of 
the SN energy heats the intercloud gas, and that the remain- 
ing ?7cid — 5% goes into the veloc ity dispersion of the cloudy 
medium (McK ee & Ostrikeilll977l) . Thus 

^icm = Hem ■ " s « • 10 51 erg * r] icm ■ p sf ■ E m , (9) 
field = Held ■ «sn • 10 51 erg w i7 cM ■ p sf ■ E m . (10) 

In the second half of these equations we have related the 
heating rates to the star formation rate; here E sa = 2.7 ■ 
10 5 pc 2 Myr~ 2 is the specific SN energy (per mass of stars 
formed for one supernova II to occur). The parameters j]\ cm , 
?7 c id and E sn are not free, but are fixed in a narrow range by 
observations and supernova models. 

(iv) SNe are the most important sources of heavy elements. 
The lack of self-consistent explosion models of massive stars 
causes uncertainties in chemical yields by factors of two and 
more. However, the detailed chemical composition is not im- 
portant for the galactic evolution. We therefore use a fidu- 
cial chemical element which traces the fraction of heavy el- 
ements produced in the hydrostatic burning phases and dur- 
in g the SN explosion of ma ssive stars. With the y i eld ta bles 
of IWooslev & Weaverl i l 19951) andlThielemann etafl dl996l) . or 
the yield approximations of Samland ( 1998) we can convert 
the abundances of the fiducial SN II element into real chemical 
abundances. This method has the advantage that it is not nec- 
essary redo the dynamical simulation when new yield tables 
become available. 



3.4.4. Intermediate mass stars 

Intermediate mass stars act as a mass storage for the ISM. The 
mass loss of these stars is significant only at the end of the stel- 
lar evolution, when these stars enter the AGB and planetary 
nebulae phase. Intermediate mass stars end as white dwarfs 
with masses between 0.5 Mq and IMq. Unlike the massive 
stars , the intermediate mass stars eject on ly weakly enriched 
gas Jvan den Hoek & Groenewegerl Il997t) and, in total, they 
return twice as much mass to the ISM as the massive stars, 
however with a much lower energy and a significant time de- 
lay. In the model we neglect the radiation of the intermediate 
mass stars, because they do not heat the ISM efficiently, but 
include the mass return and the enrichment in terms of a third 
fiducial element. 

3.4.5. Radiative cooling of the hot gas 

The heating of the ISM by SNe is mainly balanced by radia- 
tive energy losses. This cooling process determines the tem- 
perature and pressure of the hot ISM and by this has influ- 
ence on the dynamics of the ISM. We use the metallicity 
and temperature dep e ndent cool ing functions MT\ r . m ,Z\ rm ) o f 
Dal garno&McCravl ll 19721) and ISiirherland & Donital ill 99^1). 
These cooling functions provide lower limits to the real en- 
ergy loss, because they are calculated for a homogeneous gas 
in thermal equilibrium. In a real ISM with density fluctuations, 
the cooling rate can be en hanced by a factor 77 coo i = 2.3 - 10 
jMcKee & Ostrikeilfl977l) . so that 

Eicm — '/cool ' Pj cm A(7^ cm , Z[ cm ). (11) 

We use an enhancement factor of ?7 coo i = 5 to describe the ra- 
diative energy loss in the simulations. The radiative energy loss 
fij cm scales with the square of the gas density p; cm because the 
ions are excited by collisions. 



3.4.3. Supernovae type la 

The progenitors of Type la SN are believed to be close bi- 
nary stars consisting of an intermediate mass star and a white 
dwarf. In the present model only a small fraction of the 1 .5 Mq 
- 8.0 Mq stars explode as SNe la. This fraction is determined 
by the numbe r ratio of S Ne type la to type II. We take a ratio 
of 1:8.5 JSamlandlll998^ which is consistent with the observed 
SN rates in galaxies jTammann et all 1 19941) and can explain 
the observed iron abundances and the Galactic evolution of the 
[a/Fe] ratio. In our model both the energy released by a SN la 
(10 51 erg) and the mass of the white dwarf progenitor (0.6 Mq) 
are given to the hot ISM. Because the mass return of SN la is 
small compared to other stars and the total energy released is 
an order of magnitude smaller than that from SNe type II, the 
SNe la do not influence the dynamical evolution significantly. 
However, SNe type la are important for the chemical evolution, 
because they are a main source of iron-peak, r- and s-process 
elements. In the same way as for the SN II, we include a fidu- 
cial chemical element to trace the enrichment with SN la nu- 
cleosynthesis products. For the co nversion into real abun dances 
we use the chemical yield table o flNomoto et alJ (1984). 



3.4.6. Dissipation in the cloudy medium 

The cloudy medium is described as a hydrodynamic fluid char- 
acterized by its density, momentum, and kinetic energy. The 
individual clouds are assumed to be at a thermal equilibrium 
temperature. The kinetic energy density of the cloudy medium 
has source and sink terms from the heating by SNe and dissi- 
pation by cloud-cloud collisions. The description of the cloud- 
cl oud coll i sions is based on the inelastic cloud collision model 
oflL arson ( 1969). There, the energy loss of the cloudy medium 
depends on the effective cross section £ = c co n ■ nR 2 , which dif- 
fers from the geometrical cross section by the factor c co \\. For 
a medium with density p c id and velocity dispersion <x cld con- 
sisting of clouds of mass M and radius R, the kinetic energy 
density then changes according to 

8V^ 2 2 3 
£cld = -Ccoll—— 'PcldO-cld- ( 12 ) 

With the cloud mass-radius relation of Elmegreen ( 1989) we 
obtain 

field ^coll-fi^ 2 - Pepcid d3) 
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with 77 co n = 0.025 for c co u = 1 (geometrical cross section). 
Unfortunately, 77 co n is not well constrained, because magnetic 
fields, gravitational focusing, and self-gravity can influence the 
effective cloud cross section significantly. 77 co n is the most un- 
certain parameter in our ISM description. 

3.4.7. Evaporation 

The different phases of the ISM exchange mass by evaporation 
as well as condensation and cloud formation processes. These 
processes depend mainly on the sizes of the cold clouds and the 
density and temperat ure of the hot ambient gas. As in the two- 
dimensional models (Samland etall [l997). we use the model 
o flMcKee & Begelmanldl990l) and lBegelma n & McKee (1990) 
for the description of the evaporation, 



Pevap=4.35-10- 16 ^p cld -r^ 2 . 



(14) 



For 10 5 Mq clouds in the simple spherical cloud model of 
Section l3~3l we obtain an evaporation rate of 



A _ „ p- l ' 4 n T 
Hevap ~ '/evap^4 f J c\d 1 ■■ 



5/2 
icm 



(15) 



with ?7 evap = 1.0 x 10 . The value of ?7 evap depends on the 
cloud model (cloud shapes, mass spectrum). 

3.4.8. Condensation and cloud formation 

In small clouds, heat transfe r by conduction is more e ffi- 
cient than cooling by radiation IMcKee & B egelman. 1990J). In 
large clouds (cloud radius exceeds the Field-length) the cool- 
ing dominates and the clouds can gain mass by condensation. 
The timescale for condensation i s of the order of the cooling 
time of the ambient hot medium JMcKee & Begelmanl fr990). 
In addition, clouds can form when the hot gas becomes ther- 
mally unstable. This cloud formation process is important in 
high density regions (e.g., shocks) or in regions with no or low 
star formation activity. Since it is a cooling instability, this pro- 
cess also has a timescale which is approximately the cooling 
time. We use the following parameterization for the condensa- 
tion and cloud formation process: 



Pcond — ?7cond 



Picm 
^cool 



(16) 



It turns out that effectively ?] coa< i is not a free parameter. 
Reasonable values for 77 con d are in the range of 1.0 (cloud 
formation time = cooling time) to 0.3. Tests showed that if 
^cond < 0.3, the hot gas can cool to temperatures of less than 
10 4 K before condensation and cloud formation are efficient. 
We use a value of ?7 con d = 0.5, which guarantees a stable hot 
gas phase. 

3.4.9. Self-regulation 

The interaction network describes a self-regulated system in 
which moderate changes of the parameter s do n ot alter the re- 
sults significantly (see also Kopp enatall fl995). Star forma- 
tion is one example for such a self-regulation. An increase of 
the efficiency ?/ s f by a certain factor will not increase the star 



formation rate by the same amount, because the ISM will be 
heated by the stars, which in turn decreases the amount of avail- 
able cold gas to form stars. Another example is the phase tran- 
sition between hot and warm gas. A rise in the heating rate of 
the hot gas leads to an enhanced evaporation of clouds. This 
increases the density of the hot gas and enhances the cooling. 
The temperature of the gas increases only by a small amount. 
If on the other hand the heating rate decreases, condensation 
and cloud formation set in. As a result, the hot gas density and 
therefore the cooling efficiency drops and the gas temperature 
will stay nearly constant. 

Some interactions are self-regulated by dynamical pro- 
cesses. For example, a higher star formation efficiency delays 
the settling of the clouds in the galactic disk, so that the av- 
erage density of the cloudy medium is lowered which in turn 
decreases the star formation rate again. For the moment we 
neglect this (important) dynamical self-regulation, discussing 
now only the influence of the important model parameters on 
the velocity dispersion of the clouds <x c i d , the hot gas temper- 
ature T lcm , and the mass ratio between gas and clouds. We do 
this by calculating equilibrium conditions first numerically and 
then in an analytical approximation. 

The interaction network is a closed system of differential 
equations, which can be solved numerically to find equilibrium 
states for given densities of the ISM. The results are plotted in 
Fig. 0] As expected for a self-regulated system, we find that 
even density variations by five orders of magnitude have only 
small effects, changing cr cld (Ti cm ) by factors of 3 (7). 

For the analytical approach, we assume that, most of the 
ISM mass is in the cloudy medium and that the cloudy medium 
is in pressure equilibrium with the surrounding gas. Balancing 
the heating by SNe with the dissipation by cloud-cloud colli- 
sions and the cooling by radiation, and assuming that evapora- 
tion rates are equal to condensation rates, we get the following 
three relations: 

/ \ 4 / 5 

/?7 s f?7cid\ _i/io M _, 



/ \ 2 / 7 / \ 2 / 7 

TicmCC MrH (W-A 1,7 

\ ^coll / \ ?7evap / 



Picm / ?7sf ^cld 



Pcld \ 77cc.ll 



3/10 



(-) 

\ 77cool / 



1/2 



■P 



-7/20 



(18) 
(19) 



Here p is the average density of the ISM (clouds and hot gas). 
These equations confirm that there is only a weak dependence 
of cr c i d and Ti cm on the average ISM density. They are shown in 
Fig-Sfbi" comparison with the numerical results. 

In the previous discussion, we neglected metallicity effects. 
However, metallicity can be important because it increases the 
cooling efficiency. The dashed lines in Fig. |4] show the nu- 
merical equilibrium relations for an extremely metal-poor gas 
([Fe/H] = -4). The effects of low gas metallicities are similar to 
a low cooling efficiency rj mo \. In both cases, the gas temperature 
Ticm stays nearly constant, because the lower cooling efficiency 
is balanced by a higher gas density pi cm . The self-regulating 
character of the interaction network again is obvious. The equi- 
librium system is in a balance of forward and reverse reaction 
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Fig. 4. The three panels show the density dependence of the 
intercloud to cloud mass ratio, the cloud velocity dispersion 
cr c id, and the hot gas temperature in a closed-box model 
where dynamical effects are neglected. The dotted and dashed 
lines show the [Fe/H] = (solar metallicity) and the [Fe/H] = -4 
cases and the full lines show the slopes given by the analytical 
approximations of eqns. U7tlT^l 



rates. Any stress that alters one of these rates makes the system 
shift, so that the two rates eventually equalize (Le Chatelier's 
principle). 

In spite of the self-regulating character of the interaction 
network it is necessary to solve the full interaction network dur- 
ing a simulation, because processes on different timescales are 
involved that compete with the dynamical evolution. For ex- 
ample, heating processes or cloud dissipation may take longer 
than dynamical changes. Equilibrium solutions of the interac- 
tion network can be used only to study the general behaviour 
of a star-gas system where dynamical processes are slow, but 
would not be able to describe, e.g., the outflow of hot gas. 



3.4.10. Sensitivity to free parameters 

Most of the parameters in this multi-phase ISM description 
(77's) can be constrained either from theory or from observa- 
tions. The two most important parameters in the model are 
the star formation efficiency 77^, and ?7 co ii which controls the 
dissipation rate. These parameters are of special interest be- 
cause they both influence the velocity dispersion of the cloudy 
medium and thus the settling of the baryonic matter. 

In Table^we list the parameters of the interaction network 
together with their upper and lower limit values. These are dis- 
cussed in the subsections above. For rj evap and ^ co n, where there 
is little constraint, we have used a large range to illustrate that 
even in this case the resulting changes to the equilibrium are 
not large. It is obvious from the analytical model above that 
Vsf Oldd, 7coii) have influence on cr cld , r icm andp icm , while 77 ev a P 
(7/cond) can change only T Km , and by varying t] coo \ one can shift 
the gas-to-cloud mass ratio. Using the extreme values for the 
parameters given in TableQJwe obtain the uncertainties of cr c id, 
T[ cm , and pj cm (column 5, 6, and 7 of Table^. This shows that 
even large changes of the parameter values lead only to mod- 
erate variations in cr c id, T[ cm , and pj cm . From these numbers we 
conclude that <x c i d , Ti cm , and pi cm /Pcid are accurate to a factor of 
2 in our model. 



3.5. Gravitation 

Gravitation is a long range force that couples the baryonic and 
the dark matter and it is the main driver of the galactic evolu- 
tion. To determine the gravitational potential, we first map all 
stars, the different phases of the ISM and the dark matter onto 
one grid. Then, depending on the kind of grid (equidistant or 
logarithmically spaced) we solve the Poisson equation with a 
FFT (fast Fourier transformation) or a SOR (successive over- 
relaxation with Chebeyshev acceleration) method. The gravita- 
tional forces are calculated from the potential at every point. 

3.6. Initial conditions and code set-up 

The present simulations used a logarithmic 8 1 3 grid with a spa- 
tial resolution of 370 pc in the inner parts and 5 kpc in the outer 
halo. The timestep lengths, which in most cases are limited 
by the interaction processes, are of the order of 10 4 - 10 5 yr. 
The simulations start at a redshift z = 4.85 corresponding to 
t = 1.2 Gyr in a ACDM cosmology. Initially the dark halo has a 
mass of 2. 1 x 10 10 M0, distributed over a 30 kpc sphere. During 
the evolution, the halo grows by accretion to a size of 500 kpc 
with an enclosed mass of 1.8 x 1O 12 M0. The baryonic-to-dark 
matter ratio of the accreted mass is fixed at 1:5, which finally 
amounts to 3 x 10 11 Mq of baryonic matter inside the r2oo ra- 
dius. The accretion rate, taken from cosmological simulations 
dKauffmann et all Il999l see Section [3] and Fig|5} decreases 
from 32M yr" 1 at z = 4.85 (1.2 Gyr) to 8.1 Moyr 1 at z = 
(13.5 Gyr). The baryonic matter outside the r2oo radius consists 
of ionized (T > 1Q 4 K) primordial gas at the virial temperature. 
The cooling time of this gas is short compared to the collapse 
time. On this time-scale clouds form which can then dissipate 
kinetic energy and collapse inside the dark halo. Thereafter the 
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parameter 


eqn. 


range 


std. value 




T" / 'rstd 
1 lcm ' 1 icm 


Picm / Pf^ 


»7sf 


i 


(5 - 30) x 10~ 4 


5.3 x 10~ 4 


1.0-4.0 


1.0-2.7 


1.0-4.0 


77cld 


EH 


0.02-0.1 


0.05 


0.5-1.7 


0.8-1.2 


0.8-1.2 


?7cool 


EH 


2.3-10 


5 


1.0 


1.0 


1.5-0.7 


Vcon 


EH 


0.013-0.5 


0.025 


1.7-0.1 


1.2-0.4 


1.2-0.4 


^7ev;ip 


M 


(0.1 - 10) x 10- 19 


io-" 


1.0 


1.9-0.5 


1.0 


?7cond 


M 


0.3-1.0 


0.5 


1.0 


0.9-1.2 


1.0 



Table 1. Parameters of the model interaction network and the equations in which they are defined. Column 3 shows lower and 
upper limits for the parameters, column 4 the standard values used in the simulations, and columns 5, 6, and 7 give the equilibrium 
<7"cid, T lcm , and pj cm resulting with the maximum and minimum parameter values, normalized by those obtained for the standard 
model values. 



4 3 2 1 redshift o 



10 







. II II 1 1 1 


1 . 


r=r 




200 






r 'W ' : 




20 kpc u 




1 , , , 1 , , , 1 , , , 1 , , , 1 , , , - 



' I I I I I I I I I I I I I I I I I I I I I I I I I I 

2 4 6 8 10 12 14 
time [Gyr] 



Fig. 5. Baryonic mass infall rate (MIR, full line) into inner 
20 kpc and across the r2oo radius (dashed line). 

collapse is only delayed by stellar feedback processes. We have 
confirmed that starting with a cloudy gas medium leads to a 
very similar evolution. Initially, the accreted baryonic matter, 
like the dark matter, is in spherical rotation (eqn[0 with a spin 
parameter A = 0.05. It grows from inside-out and the mass 
with the highest specific angular momentum is accreted at late 
times. 

4. The formation history 

The present simulation describes the formation and subsequent 
evolution of a disk galaxy in a slowly growing dark halo. The 
halo grows by accreting dark and baryonic matter in spherical 
shells with an accretion rate that i s derived from cosmolo gical 
simulations in a ACDM universe ( Kauffm ann et all 1 19991 see 
Section [3~TT >. with a baryonic mass fraction of 1/6. The dark 
and baryonic matter have initial A = 0.05 and specific angu- 
lar m omentum distribution similar to that from Bul lock et alJ 
(200l]). Star formation, feedback, and gas physics of the bary- 
onic component are described within a two-phase ISM model 
(see Section EOl . The dynamics of the gas and stars is treated 



separately. The simulations start at z — 4.85 it = 1.2 Gyr) 
with a mixture of dark matter and gas, but no stars. At red- 
shift z = 4.85 the dark halo mass exceeds 2.0 x 10 10 M©, which 
is twice the mass resolution of the cosmological simulations on 
which our halo model is based. 

4.1. Collapse and star formation history 

The accreted primordial gas quickly reaches an equilibrium be- 
tween the cold and hot phases with most of the mass in the cold 
clouds. It can then dissipate its kinetic energy and begins to 
collapse inside the dark halo. The strong collapse takes place 
in a redshift interval between z = 1.8 (3.6 Gyr) and z = 0.8 
(6.6 Gyr). With a duration of ~ 3 Gyr, its period exceeds the 
free-fall time at the r2oo radius by a factor of 3.4. This late and 
delayed collapse, which is in contradiction to simple collapse 
scenarios, is caused by stellar feedback in conjunction with the 
initially shallow gravitational potential of the halo. 

The baryonic mass infall rate (MIR) across the r2oo radius 
and into the inner 20 kpc is shown in Figure [5] At early times 
the gas fallen through r2oo has not yet reached the inner 20 kpc 
because of the delay from feedback. With the growing mass of 
baryons the dissipation increases and the infall into the centre 
accelerates, until the combination of feedback from the grow- 
ing star formation rate and decreasing infall rate through r2oo 
reverses this trend and the MIR into 20 kpc begins to decrease 
again at z — 1.2. At late times (z < 0.3) a significant fraction of 
the infalling gas has too much angular momentum to arrive in 
the central 20 kpc, so that the rate of infall into 20 kpc becomes 
lower than that through r^oo. The mass of the galaxy increases 
rapidly and at z = 1.2 (5 Gyr) already 10 u Mq of baryonic 
matter has been accreted. At z = 0, 2/3 of the total baryonic 
mass of 3 x 10 11 Mq has collapsed into the inner 20 kpc of the 
dark matter halo and has formed the luminous galaxy, leading 
to a dark-to-baryonic matter ratio of 3%, 10% and 30% inside 
galactocentric radii of 1, 3, and lOkpc, respectively. 

The star formation rate integrated over the central 20 kpc 
is shown in Figure |6] This increases steeply to a maximum of 
~ 5OM yr _I at z = 1 (5.75 Gyr). At that time, the gas con- 
sumption by the star formation process is balanced to 80% by 
infall and to 20% by the stellar mass return. Afterwards the 
infall rate decreases to 5M©yr , at z — 0, while the stellar 
mass return stays nearly constant at a rate of lOMoyr -1 . This 
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Fig. 6. Star formation rate (SFR) in the inner 20 kpc of the dark 
halo. The star formation rate peaks slightly later than the MIR 
(Fig.0. 

explains the different shapes of the curves in Figs. [5] |6] at late 
times. 

4.2. Sequential formation of halo, bulge and disk 

Figures [3] show the spatial distributions of the gas and stars 
as a function of redshift. In general terms, star formation out of 
the dissipating cloud medium proceeds from halo to disk and 
from inside out. Thus the earliest stars at redshifts z > 2 form in 
the whole volume limited by the r2oo radius (halo formation), 
with a concentration to the central bulge region. At a redshift of 
z = 1.3 (4.7 Gyr), a (thick) disk component first appears. Later 
the star formation is concentrated to the equatorial plane, and 
at late times most of the star formation occurs in the outer disk. 

Around z = 1.2 the infall of gas can no longer compen- 
sate the gas consumption by star formation in the inner disk. 
This is because the infalling molecular cloud medium now has 
higher specific angular momentum, so that its infall timescale 
becomes longer than the central star formation timescale. As a 
result, the region of highest star formation density moves radi- 
ally outwards from the centre and a ring (Fig.0 Z = 1 panels) 
forms. This ring grows to a radius of ~ 4 kpc before, at redshift 
z = 0.85 (6.4 Gyr), the disk becomes unstable. Within 150 Myr 
the ring then fragments, and for a short time (« 150 Myr) a 
very elongated bar is formed with exponential major-axis scale- 
length of 5.9 kpc and axis ratio of 12:2:1. Fig. [5] shows that 
after this transient period the bar length decreases, and a bar- 
bulge with axial ratios ~3: 1.4:1 develops which includes the 
old bulge component formed in the early collapse. 

Already during this bar formation process the galaxy starts 
to build up the outer disk. This disk is the youngest component 
in the model galaxy, even though the oldest disk stars are as old 
as the halo stars. The disk grows from inside-out, because the 



early accreted mass has low specific angular momentum (see 
the last columns of Fig.0and Fig. [8}. In parallel, the vertical 
scale-height in a fixed radial range decreases; this is shown in 
Fig.^|for the range 3 - 10 kpc. The more pronounced settling 
of the cloudy medium to the equatorial plane at late times is due 
to (i) more efficient cooling because of higher metallicity, (ii) 
higher gas density and, thus, dissipation rate for gas near the 
angular momentum barrier in a deeper potential well. The final 
disk thickness in Fig.^|of * 500 pc is set by the resolution the 
present model. 

After the bar-bulge has formed, two trailing spiral arms ap- 
pear in the disk which are connected to the bar-bulge. At the 
beginning the two spiral arms are symmetric and of the same 
size. However, in the further evolution a persistent lopsided- 
ness of the galaxy develops. The resulting off-centre motion of 
the bulge-bar brings one arm nearer to the bulge-bar, and the 
more distant spiral arm becomes more prominent over a wider 
winding angle (Fig.0. The lopsidedness may be produced by 
the dark halo which is here described as an external spheri- 
cally symmetric potential component, even though in the inner 
galaxy the baryonic matter dominates the gravitational poten- 
tial. The effect that su ch per turbations c an produce lops ided 
galax ies is well known jLevine & Sparkeiri998HSwaters et all 

We expect this evolutionary sequence to be more or less 
typical of any dissipative collapse in a growing dark matter 
halo, unless it is interrupted by substantial mergers. The fol- 
lowing modelling uncertainties may further modify the evolu- 
tion as described above, (i) The assumed dissipation parameter 
is uncertain. If in reality dissipation is significantly more effi- 
cient than in the model, a nd in addition star s form only above a 
certain density threshold (K ennicutn. il 998). the gas would fall 
into the disk much more rapidly without significant star forma- 
tion. In this case formation of stars in the halo and thick disk 
may be substantially suppressed. Also, in this case the veloc- 
ity dispersion of the cloudy medium in the disk would decrease. 
Then the gaseous disk may fragment in local instabilities before 
a global disk instability sets in, so that more material would 
dissipate to the centre and contribute to the bulge, (ii) The pre- 
cise angular momentum distribution (as compared to the uni- 
versal average distribution assumed) will influence the detailed 
inside-out formation of the disk. For example, the formation of 
a ring in the present model would likely be suppressed by in- 
creasing the amount of low angular momentum gas falling in 
at redshift z - 1. (iii) When small-scale initial fluctuations lead 
to the formation of sub-galactic fragments that merge before 
the main collapse, the mass of the bulge component would be 
substantially greater, and its d ynamical structure d etermined by 
the merging process (Stei nmetz & Navarro! Eo02). 

4.3. Global chemical evolution 

The chemical evolution of the model proceeds in three steps. 
First, there is the collapse phase, lasting until z — 0.6. Then two 
quasi-equilibrium phases follow in which first (z = 0.6 - 0.3) 
the mass infall and later (z = 0.3 - 0.0) the stellar mass return 
determines the star formation rate (see Figs.|5]|6). These three 
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Fig. 7. Face-on surface density of the ionized gas, the cloudy medium, the OB-associations and the stars at different redshifts. 
Each column shows the evolution of one component between redshift z — 2 and z — 0. Each panel has a size of 50 x 50 kpc. 




Fig. 8. Same as Fig. [7] but for the edge-on view. 
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Fig. 9. Evolution of the bar-bulge axis ratios, z indicates the 
vertical direction and x and y the major and minor axes of the 
bar. The bar forms from a disk instability at time 6.4 Gyr and 
settles to a triaxial bulge with axis ratios ~3: 1.4:1 at late times. 
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Fig. 10. Evolution of the vertical scale height of the stellar disk 
between radial distances 3 - lOkpc. This part of the disk forms 
later than the average of all stars shown in Fig. |5] The scale 
heights in small rings are shown as dots and the radial average 
is drawn as a full line. The dashed line indicates the spatial 
resolution limit of 370 pc. 

phases are clearly visible in the chemical enrichment history of 
the cloudy medium. This is shown in Fig. where we have 
plotted the Zn/H metallicity normalized to the solar value be- 
cause Zn is believed to not be depleted onto dust grains. The 
Zn/H abundance was reconstructed from the fiducial SN I and 
SN II elements as described in Section[3] 

The initial collapse of the baryonic material produces a first 
generation of stars which already at z = 2 has synthesized 
enough metals to increase the average metallicity to ss -1. 
During the late collapse phase the star formation rate exceeds 
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Fig. 11. Average Zn/H abundance (linear scale) of the ISM and 
the stars in the inner 20 kpc of the model galaxy. 



the MIR and the metallicity increases fast. Then in the second 
phase, the star formation rate is proportional to the MIR and 
the metallicity increases only slowly. In the final phase, the star 
formation rate is again higher than the MIR and the metallicity 
increases more rapidly (Fig. II 1> . 

The chemical enrichment histories of the hot gas, cloudy 
medium, and stars (shown separately in Fig. II 11 are very dif- 
ferent from the predictions of multi-phase closed box models. 
In these models the hot ionized gas always has the highest 
metallicity, followed by the molecular cloud medium and the 
stars. This is different in our model, because stars, clouds and 
gas have different spatial distributions. Most of the ionized hot 
gas is located in the halo, while the molecular clouds and stars 
are preferentially found in the disk and bulge, where the gas 
densities and metallicities are high. After a brief initial period 
the average metallicity of the stars exceeds that of the cloudy 
medium, because the stars form from gas that is more concen- 
trated to the equatorial plane and thus more metal-rich than the 
average cold gas medium. 

An interesting feature in Fig. m is the constant [Zn/H] 
» 0.4 of the hot gas from redshift z — 1 until the present epoch. 
A large fraction of this gas occupies the galactic halo, where the 
star formation rate is low. Two processes keep the halo metal- 
licity at a constant level: Gas flows from the bulge and the disk 
transport heavy elements into the halo. There it mixes with the 
existing ISM and with infalling low metallicity gas. In this way 
the metallicity in the halo gas can remain constant over a long 
time. 

A similar enrichment can be observed in the disk. During 
the collapse the star formation and the metal production is con- 
centrated to the inner galaxy. From there, the metal-rich gas 
expands into the disk. This disk pre-enrichment causes a lack 
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of metal-poor stars ("G-dwarf problem" 1 , see Section|5|i, simi- 
lar to that observed in the Milky Way. 

We find a global age-metallicity relation and also a [O/Fe] 
to [Fe/H] relation, which are consistent with Milky Way ob- 
servations. These relations are much more influenced by the 
choice of stellar yields, resp. the star formation history, than 
by the details of the galactic model, a nd have been f ound 
in previous galactic models as well (e.g. ISteinmetz &Mulleil 
ll994tlTimmes et allll995USamlandlll998l) . However, the pre- 
enrichment of the galactic disk, explaining the so-called "G- 
dwarf problem" and the sub-solar abundances in the outer 
halo gas are a result of dynamical processes and are therefore 
strongly influenced by the galactic model. The same is true for 
relations connecting the kinematics and chemical composition 
of stars; see Section l5"71 

In Fig. we compa re our results with the [Zn/H] data 
for d amped Lya- systems JPettini et allll999t lp"rochaska et all 
12001 1) , using [Zn/H] because [Fe/H] may be depleted by dust. 
Damped Lye systems are an important test for galactic models. 
If they are associated with forming disk galaxies, they can be 
used to measure the metallicity of the ISM of young, gas rich 
galaxies at high redshifts. The Lya systems on average show 
only a mild evolution in [Zn/H] for z ~ 2-3. Even though 
we consider only a single model with its specific evolutionary 
history, the predicted metallicity range as a function of redshift 
is consistent with the data. From the model we thus infer that 
the time needed to increase the disk metallicity to the typical 
values observed in Lya systems is ~ 1 Gyr. The shaded area 
in Fig.ll2lshows the range of metallicities from different loca- 
tions in the model at redshifts z > 1.6. Fig. El also shows the 
early [Zn/H] evolution of the hot ionized gas (dashed line). This 
should be representative for the average metallicity of metal- 
enriched clouds forming in the low density halo ISM, and thus 
for absorption of background quasars along halo lines-of-sight. 

4.4. The ISMatz = 

Atz = O,5.6xlO lo M (15.7%) of the initial baryonic mass of 
3 x 10 11 Mq is still in the form of gas. The distribution of this 
late ISM is spatially more extended than that of the stars and 
the mass ratio of ISM to stars increases with distance from the 
galactic centre. Inside lOkpc, only about 2% of the baryonic 
mass is still in the ISM, but inside 60 kpc the fraction is already 
43%. 75% of the ISM mass is in the cold and warm phases 
(molecular and atomic gas), but a significant fraction of 25% is 
hot and ionized. This hot gas extends far out into the halo and 
fills most of the halo volume (see Fig. [3}. As discussed in the 
previous subsection, the hot gas is the chemically least enriched 
component at late times. 

5. Formation of galactic components 

The model presented in this paper describes the formation of a 
massive disk galaxy with a total luminous mass of 2 x 10" Mq 

1 The term "G-dwarf problem" describes the fact that the fraction 
of low-metallicity stars in the Milky Way disk is substantially lower 
than predicted by simple 1-zone models. 
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Fig. 12. Redshift evolution of the average zinc abundance 
[Zn/H] of the cloudy medium (full line) inside a galactocen- 
tric radius of 20 kpc. The shaded area shows the range of 
[Zn/H] metallicities in the model ove r this region. The p oints 
are data for damped Ly g systems from lPettini et alJ (Tl999) and 
Prochask a et alJ (12001 ). plotted as squares and rhombi, respec- 
tively. The general trend of the model curve is consistent with 
the observed damped Lya abundances. The dashed line indi- 
cates the metallicity evolution of the hot gas component, which 
fills most of the galactic halo. 

inside 20 kpc and a disk rotation velocity of 270 km s _I at a ra- 
dius 10 kpc, about 2.5 times more massive than the Milky Way. 
The galaxy model forms from inside-out radially and from halo 
to disk vertically. In Fig. |6] one can clearly see halo, bulge, 
and disk components, which as described in Section 0] form 
sequentially in this order. The model contains the full spa- 
tially resolved dynamical and chemical information for all of 
these components. With this information it is possible to study 
the detailed consequences for the stellar kinematics and abun- 
dances, that are implied by the assumptions made for the galaxy 
formation scenario and the modelling of the various physical 
processes involved, and to compare to a variety of observa- 
tional data. 

Figure [O] shows average properties of the model stars at 
redshift z = 0: age, metallicity, [O/Fe], and rotational velocity. 
These vary continuously as a function of position in the merid- 
ional plane (R, Z), in a way that reflects the interplay between 
star formation and enrichment on the one hand, and the verti- 
cally top-down and radially inside-out dynamical formation on 
the other hand. Thus, the oldest stars can be found at low radial 
distances in the halo and the youngest stellar populations are 
located in the outer galactic disk. Because mean age and mean 
metallicity vary as a function of position, global mean values 
for age or metallicity are not sufficient to characterize the stel- 
lar population of a galaxy. Notice also that stars at a given mean 
age or metallicity can be found in a variety of positions. 
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Fig. 13. The panels show the mean stellar age, [Fe/H], [O/Fe], and mean rotation velocity V of the model stars at redshift z — 0, 
in the meridional plane as a function of radial (R) and vertical distance (Z) from the galaxy's centre. The mean stellar ages are 
given in Gyr, and the units of the velocities are pc Myr 1 = 0.98 km s _1 . 



In Fig. [2] we plot the complete distribution of model stars 
at redshift z — in the plane of rotational velocity v lot ver- 
sus metallicity [Fe/H]. This shows the expected evolution from 
non-rotating, low-metallicity stars formed early in the collapse 
to a rotating disk-like, metal-rich population formed at late 
times. However, this evolution is by no means described by a 
single dependence of v rot on [Fe/H]; rather there is large scat- 
ter in v rot at given metallicity, especially at low [Fe/H] values. 
In addition, there are also at least two concentrations of stars 
visible, at [Fe/H] a -0.4 and [Fe/H] ^0.1. 

The top panel of Fig. 1141 shows a [Fe/H] histogram of all 
stars, projecting along the v rot axis. Based on this histogram, 



we will now study the stellar populations in the different parts 
of Fig. ^] separated by the dashed lines. In first order these 
might be described as extreme halo, inner halo, metal-weak 
thick disk, thick disk, thin disk, and central bulge components, 
somewhat reminiscent of the similar components known in the 
Milky Way. However, we do not intend to imply a detailed cor- 
respondence with the Milky Way components; there are signif- 
icant differences with the present model. Recall that the dark 
and baryonic masses of this model are larger than in the Milky 
Way, and especially the dissipation parameter of the model may 
not be typical for the Milky Way. Nonetheless, studying these 
components is useful to understand the course of the dynamical 
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[Fe/H] 

Fig. 14. Lower panel: Distribution of model stars in the 
metallicity-rotation velocity plane. Rotation velocities are per- 
pendicular to the total angular momentum vector and are in 
pcMyr . Upper panel: total metallicity distribution of all 
model stars, projecting along the rotation velocity axis. The 
dashed lines separate the subpopulations discussed in the text. 

evolution and chemical enrichment in dissipative, star-forming 
collapse, and we will for definiteness use the Milky Way ter- 
minology in the following discussion of the properties of these 
populations. 

5.1. [Fe/H]< -1.9 ("Extreme halo") 

This component comprises all stars with [Fe/H] < —1.9; these 
stars form before the collapse and simultaneous spin-up of the 
gas goes under way. As the left row of Figure[2]shows, these 
stars form mostly in the first two Gyr after the start of the 
simulation (at time 1 .2 Gyr), but there is a small ongoing star 
formation in the halo even to late times from infalling low- 
metallicity gas (remember that we have not imposed a density 
threshold for star formation in this model). These stars have 
[O/Fe] near 0.4 typical for enrichment by SNe II. Their spatial 
distribution is nearly spherical, they rotate slowly if at all, and 
they have an asymmetric, broad distribution of eccentricities bi- 
ased towards radial orbits. Here we have defined eccentricity as 
e = (r max - r min )/(r max + r min ) where r max and r min are the maxi- 
mum and minimal three-dimensional radii of the orbit. The rel- 
ative lack of low-eccentricity orbits at these metallicities (see 
Section l5^7t may be due to our neglec t of small-scale structur e 
in the initial conditions; compare with Bekki & Chiba (2001). 

5.2. -1.9 <[Fe/H]< -0.85 ("Inner halo") 

Stars in this component have a metallicity range - 1 .9 < [Fe/H] 
< -0.85, covering the spin-up phase until the first minimum in 
the [Fe/H] histogram in Fig.^] These stars form a more flat- 



tened spheroidal distribution with a flat density profile. They 
rotate with ~ 70 km s _1 , and the eccentricity distribution is less 
radially biased than for the extreme halo component. This com- 
ponent forms partly during the phase of rapid increase in the 
global star formation rate (Fig.|6jl, but there is also a younger 
part forming until late times. Here the selection in metallicity is 
important: the older part of these stars forms at relatively small 
radii when the metallicity there reaches the selected range, 
while the younger component forms at progressively larger 
radii from infalling low-density gas mixed with outflowing en- 
riched material. This younger component may not form in real 
galaxies if stars can only for m above a density threshold, as 
inferred from observations by iKennicuttl Jl998l) . Further sim- 
ulations also showed that a larger dissipation parameter (see 
Section |3. 4. 61 leads to a faster settling of the clouds into the 
disk and this suppresses the star formation in the halo at late 
times. 

5.3. -0.85 <[Fe/H]< -0.6 ("Metal-weak thick disk") 

These stars with -0.85 <[Fe/H]< -0.6 are intermediate in their 
rotational properties between a halo and a disk population (see 
Fig. 114-b . As Fig.ll5lshows their spatial distribution is already 
strongly flattened. They have a sizeable rotation with peak of 
the distribution at v rot - 150kms~', around 50-60% of the cir- 
cular velocity. The eccentricity distribution is broad and cen- 
tred around 0.5. The oldest stars of this component form at 
around the peak of the global star formation rate when the gas 
still has [O/Fe] ^ 0.3. A slightly larger fraction of stars in this 
range of [Fe/H] forms continuously until late times, at progres- 
sively larger radii, again favoured by the lack of star formation 
threshold. These stars are already enriched by SN la, as visible 
in their [O/Fe] distribution. 

5.4. -0.6 < [Fe/H] < -0.15 ("Thick disk") 

This component with -0.6 <[Fe/H]< -0.15 is one of the dis- 
tinct components in Fig. El an d is clearly dominated by a disk 
component in its rotational properties and eccentricity distribu- 
tion. However, Fig.^^shows that it also contains intermediate 
metallicity retrograde stars, showing that this metallicity selec- 
tion also includes part of the bulge. These bulge stars cause the 
peak in the age distribution at ~ 5 Gyr, with high [O/Fe]. The 
disk part has a relatively narrow distribution in [O/Fe] centred 
around ~ 0.13. 

5.5. -0.15 <[Fe/H]< 0.17 ("Thin disk") 

The thin disk component, selected in the interval -0.15 < 
[Fe/H] < 0.17, is characterized by near-circular orbits and a 
high rotation velocity, about 270 km s _1 . This component forms 
with an approximately constant star formation rate and has 
mostly near-solar [O/Fe]. As in the previous thick disk compo- 
nent, a contribution from bulge stars is included by the metal- 
licity selection. The thin disk shows the bar and spiral arms 
clearest; see the top panels in Fig. [21 In time, it forms from 
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inside out due to the accretion of material with progressively 
larger angular momentum. 

The thin disk shows only a shallow radial metallicity gradi- 
ent of 0.02 dex/kpc, which is due to efficient mixing of the disk 
ISM by the galactic bar. Beyond 16kpc, where the mixing is 
not very efficient, a steeper metallicity gradi ent of .07 dex/kpc 
is established. See also iMartinet & Friedlil dl997l) who found 
strong metallicity gradients only in non-barred galaxies. The 
metal production by the bulge and the mixing by the bar lead 
to a pronounced pre-enrichment of the disk, which is visible in 
Fig.[H)]by the lack of metal-poor stars. 

5.6. [Fe/H]> 0.17 ("Inner bulge") 

The final metallicity component with [Fe/H] > 0.17 consists 
mainly of stars in the central bulge. Notice that these most 
metal-rich stars do not form preferentially at late times; their 
age distribution is similar to that of the thick disk component. 
As the lower panels of Fig.^]show, the metal rich inner bulge 
forms from z - 1 in this model, rotates relatively slowly and 
has a broad eccentricity distribution. 

The total bulge population includes this inner component 
and a second, more metal-poor component included in the 
metallicity range of the "thick disk". Two sub populations in 
the bu lge is consistent with the observations of Pru gniel et alJ 
who find that galactic bulges in general consist of two 
different stellar components, an early collapse population and 
a population that formed later out of accreted disk mass. 

5.7. Kinematics-abundance relations 

The previous subsections have shown that while it is useful 
to select stellar components by metallicity, such a selection 
is not ideal because it lacks spatial and kinematical informa- 
tion. Vice-versa, a purely kinematical selection would mix stars 
from halo and bulge, or metal-poor and metal-rich disk. A more 
physical separation would involve spatial, kinematic, and abun- 
dance information. This is possible with models such as this, 
because the complete population of model stars is available, 
but will be deferred to a later paper. Here, we will only show 
two relations between kinematic and abundance parameters for 
all stars in the model, which further illustrate some of its prop- 
erties. 

Figure^]shows the mean stellar rotation velocity as a func- 
tion of metallicity [Fe/H] for all model stars. For [Fe/H] < -2 
v rot is independent of metallicity, at about 50kms~'. It then 
rises roughly linearly with [Fe/H] until the disk rotation veloc- 
ity is reached at [Fe/H] =s 0.4. As is apparent from Fig. [H] 
there is a large dispersion around the mean value except in the 
disk components. Figure^Dshows the mean orbital eccentric- 
ity and a sample of individual stellar eccentricities (defined in 
Section l5~Tt as a function of metallicity. It is clear that at any 
metallicity, there is a large scatter in the orbital shapes. One of 
the causes of this is the very substantial deepening of the gravi- 
tational potential during the prolonged infall, which also affect 
orbits of stars formed early in the collapse. A second effect is 
that the turbulent velocities in the cloud medium are substantial 



during the collapse and during periods of large star formation 
rate, so that the stars made from these clouds inherit significant 
initial velocities. Thus all eccentricity distributions shown in 
Fig. ll5l are fairly broad. When the various distributions shown 
in that figure are added with the appropriate mass weighting, a 
broad distribution results. 

Both the dependence of v rot on metallicity in Fig. 1171 and 
that of eccentricity on [Fe/H] in Fig.^Jare simi lar to those ob - 
served in the Milky Way near the Sun JChiba & Beers! EoOO). 
even though our model was not made to match the Milky Way 
and does not in detail describe the Milky Way. The main con- 
clusion we draw from diagrams like Figs.EHDDis that realistic 
dissipative collapses are much more complica ted than the sim- 
ple model put forward bvlE ggen et alJ II 19621) forty years ago. 
Therefore, inferring a merger history from the lack of rotation 
gradient with metallicity and a broad eccentricity distribution 
in the Milky Way halo is premature, irrespective of the fact 
that a fraction of the h alo is likely to have been formed out 
of merging fra gments dlhntn at nil 12001 1 iHelmi etall [l999; 
IChiba & BeerslUoOOh . The important question, about the rel- 
ative importance of the disruption of small scale stellar frag- 
ments, merging of gaseous fragments, and star formation in 
smooth dissipative accretion, all within their dark matter halos, 
will need more detailed modelling and more data to answer. 
All of these processes are expected in current CDM galaxy for- 
mation theories. The present simulations show that feedback is 
important for the final properties of the stars formed. 

6. Conclusions 

We have presented a model for the dissipative formation of a 
large disk galaxy, inside a growing dark matter halo accord- 
ing to the currently favoured ACDM cosmology with Qa=0.7, 
fio=0.3, ho=0.7. The dark halo mass accretion history is taken 
from the large-scale str ucture formation simulations of the 
GIF- VIRGO consortium JKauffmann et allll999t) . We use a to- 
tal dark matter mass of 1.5 x 1O I2 M0, a total baryonic mass 
of 3 x 10 u Mq, a spin parameter A = 0.05, and an angu- 
lar momentum pro file similar to the universal profile found by 
lBullocketaild200lh . 

We model the dynamics of stars and of two phases of 
the interstellar medium, including a phenomenological descrip- 
tion of star formation, the most important interaction processes 
between stars and the ISM, and the feedback from massive 
stars and SNe. The dynamical evolution of the stars is treated 
with a particle-mesh code. For the dynamical description of 
the cold molecular and hot ionized gas components we use a 
three-dimensional hydrodynamical grid code. Linked to the N- 
body and hydrodynamical code is a Poisson solver for the self- 
gravity and an interaction network, including a description of 
star formation, the heating and enrichment of type la and II 
SNe, mass return by intermediate mass stars, the heating and 
cooling by radiation, dissipation in the cloudy medium, for- 
mation and evaporation of cold clouds in the hot intercloud 
medium. We discuss the range of uncertainty in the parame- 
ters that enter this description. Because of the self-regulation of 
these processes that is quickly established, even large changes 
in these parameters only lead to moderate changes, about a fac- 



Markus Samland and Ortwin E. Gerhard: Formation of a Disk Galaxy 



19 







# 


tu 





m 


- 


— 




f 



3 4 6 6 10 12 3 4 6 8 10 13 2 4 6 S 10 12 2 4 6 S 10 12 
time [Gyr] time [Gyr] time [Gyr] time [Gyr] 



0.0 0.2 0.4 
[O/Fe] 



0, 0.2 0.4 
[QTe] 



0.0 0.2 0.4 
[O/Fe] 



0.0 0.2 0.4 
[CWe] 




L 




0.0 0.2 0.4 0.6 0.S 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 
eccentricity eccentricity eccentricity eccentricity 



A 



r 



-200 200 400 
v [pc/Myr] 



-200 200 400 
v^^ [pc/Myr] 



-200 200 400 



-200 200 400 
v [pc/Myr] 



Fig. 15. Properties of model subcomponent selected by metallicity according to the dashed lines in Fig. El Columns from left 
to right show all stars with [Fe/H] < -1.9 ("Extreme halo"), -1.9 < [Fe/H] < -0.85 ("Inner halo"), -0.85 < [Fe/H] < -0.6 
("Metal-weak thick disk"), -0.6 < [Fe/H] < -0.15 ("Thick disk"), -0.15 < [Fe/H] < 0.17 ("Thin disk"), and [Fe/H] > 0.17 
("Inner bulge"). For each of these components, the panels from top to bottom show the face-on projection onto the disk plane, 
edge-on projection, distribution of formation times, [O/Fe] distribution, eccentricity distribution, and distribution of rotation 
velocities. The top two panels in each row show an area of 40 kpc x 40 kpc. Star formation in the model starts at time 1 .2 Gyr. 
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tor of 2, in the cloud velocity dispersion, hot gas temperature 
and ratio of hot to cold gas density. 

The simulation starts at a redshift of z — 4.85 with a dark 
halo of 2.1 x 10 i() Mq. Inside the growing dark halo a disk 
galaxy forms, whose star formation rate reaches a maximum 
of 5OM yr _1 at redshift z ~ 1. Our main results can be sum- 
marized as follows: 

1. The formation and evolution of a galaxy in a dissipa- 
tive collapse scenario is much more complex than predicted by 
simple models. The energy release by SNe and massive stars 
(feedback) prevents the protogalaxy from a rapid collapse and 
delays the peak in the star formation to a redshift of z ~ 1 . The 
prolonged formation time causes fundamental changes in the 
galactic shape, the kinematics of the stars, and the distribution 
of the heavy elements. 

2. The galaxy forms radially from inside-out and vertically 
from halo-to-disk. The dynamics of the collapse strongly in- 
fluences the star formation and chemical enrichment history, 
by modifying the gas density, cooling time, and thus star for- 
mation rate. The on average youngest and oldest stellar pop- 
ulations are found in the outer disk and in the halo near the 
galactic rotation axis, respectively. As a function of metallicity, 
we have described a sequence of populations, reminiscent of 
the extreme halo, inner halo, metal-poor thick disk, thick disk, 
thin disk and inner bulge in the Milky Way. The first galactic 
component that forms is the halo, followed by the bulge and 
the disk-halo transition region, and as the last component the 
disk forms. 

3. An interesting feature of this model is the formation of a 
gaseous ring at a redshift of z — 1 which collapses and forms a 
galactic bar. The bar induces numerous changes of the galaxy's 
properties. It alters the density structure of the disk, induces 
spiral arms, enhances the mixing of ISM, and flattens metal- 
licity gradients in the disk. It also channels gas into the centre, 
leading to prolonged star formation in the bulge region. 

4. The bulge thus contains at least two stellar populations: 
An old population that formed during the proto-galactic col- 
lapse and a younger bar population. The distinguishing feature 
between these populations is the [a/Fe] ratio, with the bar pop- 
ulation reaching [a/Fe] < 0. 

5. The disk is the youngest galactic component, charac- 
terized by an approximately uniform star formation rate at 
late times. Its metallicity is approximately solar and its stel- 
lar metallicity distribution shows a pronounced lack of low- 
metallicity stars compared to simple closed-box models, due to 
pre-enrichment of the disk ISM. 

6. Early in the collapse an approximately spherical halo 
forms. As the collapse proceeds, the dissipation leads to spin- 
up of the gas, and newly formed stars acquire increasing ro- 
tation speeds, as expected. However, the distribution of orbital 
eccentricities as a function of metallicity for halo stars has large 
scatter. As a result, this distribution and the mean rotation rate 
as a function of metallicity are not very different from those 
observed in the solar neighbourhood. This shows that early ho- 
mogeneous collapse models are oversimplified, and that con- 
clusions based on these models about the degree of merging 
involved in the formation of the Milky Way halo can be mis- 
leading. 



7. The metal enrichment history in this model is broadly 
consistent with the evolution of [Zn/H] metallicity in damped 
Lya systems. 

8. The most metal-rich stars form approximately 1 Gyr after 
the peak of the star formation, when the SN la rate is at its 
maximum. These stars have the lowest [a/Fe] ratios and, on 
average, an age of ~ 5 Gyr in this model. The infall of gas and 
the mass return from old stars decreases the average metallicity 
in the inner galaxy at later times. 

This chemo-dynamical model provides kinematics and 
metallicities of individual stars, but also can be used to ob- 
tain colours, metallicities and velocity dispersions of integrated 
stellar populations. It can therefore be used as a tool to under- 
stand observations of distant spiral galaxies, and connect them 
with stellar data in the Milky Way. In addition, it provides gas 
phase metallicities and temperat ures, and star format ion rates, 
as a function of time or r edshift. Western et al. (2002) have al- 
ready used this information to predict the colour evolution of 
large spiral galaxies, and have compared with bulge colours in 
the Hubble Deep Field. Further "observing" of such models 
and comparing with distant galaxy observations may lead to 
valuable insights on the galaxy formation process. 

A table of the positions, velocities, and metallicities for all 
stars in the model can be obtained by request from the authors. 
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Fig. 15. (continued) 
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Fig. 16. Distribution function of oxygen abundances of G- 
dwarfs, selected by mass, in the model galaxy at redshift z — 0. 
All G-dwarfs in the equatorial plane with galactocentric dis- 
tances between 9.5 kpc and 10.5 kpc are plotted. The dotted 
curve shows predictions of an instantaneous simple model with 
an effective yield Y e g - Z . 



Fig. 18. Sample of stellar orbital eccentricities and mean eccen- 
tricity of all stars in the model, as a function of metallicity. 




Fig. 17. Mean stellar rotation velocity as a function of [Fe/H] 
for all model stars. 



